The role of torsional Alfven waves in coronal heating 
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ABSTRACT 

In the context of coronal heating, among the zoo of MHD waves that exist in 
the solar atmosphere, Alfven waves receive special attention. Indeed, these waves 
constitute an attractive heating agent due to their ability to carry over the many 
different layers of the solar atmosphere sufficient energy to heat and maintain a 
corona. However, due to their incompressible nature these waves need a mecha- 
nism such as mode conversion (leading to shock heating), phase mixing, resonant 
absorption or turbulent cascade in order to heat the plasma. Furthermore, their 
incompressibility makes their detection in the solar atmosphere very difficult. 
New observations with polarimetric, spectroscopic and imaging instruments such 
as those on board of the japanese satellite Hinode, or the SST or CoMP, are 
bringing strong evidence for the existence of energetic Alfven waves in the solar 
corona. In order to assess the role of Alfven waves in coronal heating, in this work 
we model a magnetic flux tube being subject to Alfven wave heating through the 
mode conversion mechanism. Using a 1.5-dimensional MHD code we carry out 
a parameter survey varying the magnetic flux tube geometry (length and expan- 
sion), the photospheric magnetic field, the photospheric velocity amplitudes and 
the nature of the waves (monochromatic or white noise spectrum). The regimes 
under which Alfven wave heating produces hot and stable coronae is found to be 
rather narrow. Independently of the photospheric wave amplitude and magnetic 
field a corona can be produced and maintained only for long (> 80 Mm) and 
thick (area ratio between photosphere and corona > 500) loops. Above a critical 
value of the photospheric velocity amplitude (generally a few km s _1 ) the corona 
can no longer be maintained over extended periods of time and collapses due to 
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the large momentum of the waves. These results establish several constraints on 
Alfven wave heating as a coronal heating mechanism, especially for active region 
loops. 

Subject headings: Sun: corona - Sun: flares - MHD - waves 



Introduction 



New observations from polarimetric, spectroscopic and imaging instruments are re- 
vealing a corona permeated with waves. Compressive modes such as the slow and fast 
magnetoh ydrodynamic (MHD) modes have a long chapter in observational history (see re 
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views by iNakariakov fc Verwichtd 120051 ; iBanerjee et al.l 120071 : iRuderman fc Erdelyil 12009 



Only recently, the development of high resolution instruments 
has brought strong evidence for the existence of the third MHD mode, the Alfven mode, in 
the solar atmosphere. In the majority of the observational reports of Alfven waves ambi- 
guity exists and the wa ves can be interpreted as trapped modes, fast and slow kink waves. 
De Pontieu et al.l (120071 ) reported transversal displacements of spicules from observations in 
the Ca II H-line (3968 A) with a broadband filter of Hinode /SOT. The wavelengths of these 
chromospheric agents were estimated to be longer than 20000 km, with periods between 
100 and 500 s and speeds of at least 50-200 km s _1 . In the absence of a stable waveg- 
uide they interpreted the swaying of t he spicules as a result o f upward propagating Alfven 
waves. However, it was pointed out by lErdelyi fc Fedunl (120071 ) that these waves were likely 
to be kink oscillations, due to the displacement of the axis of symmetry o f the flux tube 



caused by the later. Using the Coronal Multi-Channel Polarimeter (CoMP) iTomczyk et al. 



( 120071 ) analyzed properties of infrared coronal emission line FeXIII (1074.7 nm) across a 
large field of view with short integration times. Waves propagating along magnetic field 
lines with ubiquitous quasi-periodic fluctuations in velocity with periods between 200 and 
400 s (power peak at 5-min) and negligible intensity variations were reported. Wavelengths 
and phase speeds were estimated to be higher to 250 Mm and 1 Mm s _1 respectively. In 
this work it was suggested that these waves were Alfven waves probably being generated in 
the chromospheric network from mode conversion of p-modes propagating from the photo- 
sphere (hence explainin g the peak in the power spec t rum o f the waves) . On the other hand 
Erdelyi &: Fedunl ( 20071 ) and lVan Doorsselaere et al.l ( 20081 ) argued that an interpretation in 
terms of fast kink waves was more appropriate due mainly to the observed collective behavior 
which should be absent in the case of Alfven waves. Being compressible waves it was argued 
that kink waves would appear however incompressible in the corona for an instrument such 
as CoMP due to the very small variation in intensity they produce. 
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The difficulty in detecting the Alfven wave is due in part to its incompressible nature. 
Being only a transverse disturbance in the magnetic field, it is practically invisible to imaging 
instru ments, unless t he structure of propagation is displaced periodically from the line of 
sight (jWilliamsl 120041 ) . Also, due to the large wavelengths (a few megameters) and short 
periods (a few minutes) that may be involved in the corona, a large field of view, short 
integration time and proper resolution are needed for their detection in the corona. With 
po larimeters and spectrogr aphs however Alfven waves are more easily detected, as suggested 
by lErdelyi fc Fedunl (120071 ). These waves propagate as torsional disturba nces of the magneti c 
field, causing a periodic and spatially dependent spectral line broadening (IZaqarashvilill2003l ) . 
Up to dat e , the only observational report on torsional Alfven waves is the recent work by 
Jess et al.l ( 120091 ). in which observations in Ha of a bright point using the Solar Optical 
Universal Polarimeter (SOUP) of the SST were reported. No variation in intensity nor line- 
of-sight velocity was detected and a periodic spectral broadening of the Ha line was found, 
leaving Alfven waves as the only possible interpretation for their results. The waves were 
reported to have periods in the 100-500 s interval (with maximum power in the 400-500 
s range). By comparing Ha continuum with Ha core images they reported a canopy-like 
structure with magnetic fields expanding ~ 1300 km in a vertical distance of ~ 1000 km. 
This large expansion is crucial for Alfven wave heating theories. For instance, assuming an 
Alven speed of 10 km s _1 in the upper photosphere/ lower chromosphere and a period of 
100 s (a value in the lower range of periods reported in the works cited above) we obtain 
a wavelength of 1000 km for Alfven waves propagating from the photosphere, matching the 
distance of large expansion of the flux tubes. Moreover, the height in the atmosphere where 
this takes place is the region of transition from high beta to low beta plasma. We thus have 
ideal conditions for mode conversion. Trapped modes will thus exchange their energy and a 
wide range of waves with different characteristics are likely to issue from this region. 

The large emphasis that has been put on the search for Alfven waves in the solar corona 
is largely due to their connection with coronal heating. Theoretically, they can be easily 
generated in the photosphere by the con stant turbulent conyective motions, which inputs 
large amounts of energy into the waves (jMuller et al.l 1 19941 ; IChoudhuri et al.lll993l ). Hav- 
ing magnetic tension as its restoring force the Alfven wave is less affected by the large 
transition region gradients with respect to other modes. Also, when traveling along thin 
magnetic flux tubes they are cut-off free since they are not coupled to gravity (Musielak et 
al. 2007)0. Alfven waves generated in the photosphere are thus able to carry sufficient en- 



Verth et al. ( 2009. in preparation! ) have pointed out however that the assertion made by Musielak et al 



([20071 ) is valid only when the temperature in the flux tube does not differ from that of the external plasma. 
When this is not the case a cut-off frequency is introduced. 
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ergy into the corona to compensate the losses due to radiation and conduction, and, if given 



a suitable dissipation mechanism, heat the p 
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Another possible generation mechanism for Alfven waves is through magnetic recon- 
nection. The amount of energy imparted to these waves during the reconnection process 
may depend on the location in the atmosphere of the event and is a subject of controversy. 
Parker! (Il99ll ) suggested a model in which 20 % of the energy released by r econn ection events 
in the solar corona is transfered as a form of Alfven wave. lYokoyamal (119981 ) studied the 
problem simulating reconnection in the corona, and found that less than 10 % of the total 
released energy goes into Alfven waves. This result i s simil ar to the 2-D simulation results 
of photospheric reconnection by iTakeuchi fc Shibatal (120011 ). in which it is shown that the 
energy flux carried by the slow magnetoacoustic waves is one order of magnitude higher 
that the energy flux carried by Alfven waves. On the other hand, recent simulations by 
Kigure et al. (private communication) show that the fraction of Alfven wave energy flux 
in the total released magnetic energy during reconnection in low (3 plasmas may be signif- 
icant (more than 50 %). Since the observed ubiquitous intensit y bursts (nano flares) are 
thought to play an important role in the heating of the corona (jHudsonl Il99ll ) and since 
they are generally assumed to be a signature of magnetic reconnection it is then crucial to 
determine the energy going into the Alfven waves during the reconnection process. More- 
over, it is possible that the observed na n oflare s themselves are a consequence of Alfven wave 
heating as porposed b yMoriyasu et al.l (120041 ) . In that work the resulting intensity bursts 
producing the nanoflares are created by Alfven waves by first converting to longitudinal 
mo des which steepen in to shocks and heat the p lasma. This model w as further developed 
by lAntolin et al.l (120081 ) (hereafter, Paper 1) and lAntolin et al.l ( 20091 ). where it was shown 
that the frequency of the resulting heating events from Alfven waves followed a power law 
distribution. It was further shown that Alfven wave heating could be differentiated from 
nanoflare-reco n nection heating during o bservations through a series of signatures (see also 
Taroyanl 120081 ; iTaroyan fc Erdelyil 120091 ) . For instance, Alfven waves lead to a dynamic, 
uniformly heated corona with steep power law indexes (issuing from statistics of heating 
events) while nanoflare-reconnection heating leads to lower dynamics (unless catastrophic 
cooling takes place in the case of footpoint concentrated heating) and shallow power laws. 

The main problem faced by Alfven wave heating is to find a suitable dissipation mech- 
anism. Being an incompressible wave it must rely on a mechanism by which to convert 
the magnetic energy into heat . Seyeral dissipation mech anisms have been proposed, such 



as parametric decay (I Goldst einl 1 1 9 78l ; iTerasawa et al.lll986l ). mode conversion (IHollweg et al. 
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Erdelyi &: Goossensl 1995 ). The main difficulty lies in dissipating sufficient amounts 
of energy in the correct time a nd space scales. F o r more discussion regar ding t his issue the 
reader can consult for instance iKlimchukl (120061 ); lErdelyi fc Ballail (120071 ) and lAschwanden 

(SoJ. 



Due to the observed increasing importance of Alfven waves in the solar atmosphere 
in this work we address the subject of coronal heating by Alfven waves in which mode 
conversion and parametric decay are taken as dissipation mechanisms. We concentrate on the 
heating of closed magnetic structures which populate the solar atmosphere, such as coronal 
loops, and analyze the efficiency of this heating model by carrying out a parametric space 
survey. We take a 1.5- dimensional model and we consider different photospheric drivers 
(a random driver creating a white noise spectrum, or a monochromatic driver in which 
several periods are tested), vary the photospheric magnetic field, the loop expansion from 
the photosphere to the corona and the loop length, and determine the regimes for which 
Alfven wave heating plays an important role in coronal heating. The paper is organized 
as follows. In §2] we set up the loop geometries and define the equations. As in Antolin 
et al. (2008) the loops are modeled with a 1.5-dimensional MHD code including thermal 
conduction and radiative cooling. In $3] we analyze the effect of varying the parameters on 
the thermodynamic structure of the loop. Further discussion and conclusions are presented 
in M 



2. Alfven wave model 

2.1. Parameter space 

The parametric survey in this work involves parameters changing the geometry of the 
loop, such as the loop's length and its expansion from the photosphere to the corona. We 
consider lengths of 100, 80 and 60 Mm, leading to lengths of the corona of roughly 70, 50 
and 20 Mm respectively, which are typical lengths for coronal loops. 

Due to the passage from gas pressure-dominated (high (5) to magnetic pressure-dominated 
(low (3) plasmas as we move to higher layers from the photosphere, magnetic flux tubes are 
expected to expand in the photosphere and chromosphere, forming t he so-called "magnetic 



canop y" , and have a roughly constant cross-section area in the corona (lAschwanden fc Nightingale 



20051 ). The loop expansion is a parameter influencing considerably the behavior of the waves 



as they propagate through the solar atmosphere. Apart from influencing the density strati- 
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fication, it also introduces steep gradients in the Alfven velocity due to the gradient in the 
magnetic field. The latter will in turn introduce different cut-offs for the waves propagating 
along the structure. Also, refraction of the fast MHD mode can easily occur, as fast MHD 
waves will always tend to follow the regions of higher (3 plasma. Hence, expansion can de- 
termine not only the properties but also the nature of the waves that propagate along the 
magnetic field. Here we consider different expansions of the loop displaying a cross-sectional 
area ratio between the photosphere and the corona of 1000, 600 and 500. 

Apart from the previous geometrical parameters, we also consider different photospheric 
magnetic field strengths. Observations have shown the existence of kilo-gauss photospheric 
flux tubes all over the solar surface. Valu es have been reported ranging from below 1 kG to 



above 2 kG with a mean around 1.5 kG (jSolankil Il993l ) . Here we survey different values for 



the photospheric magnetic field at the foopoint of the loop, namely, 2.3, 1.5 and 1 kG, and 
investigate the influence on the thermodynamics of the loop. 

Another parameter we consider concerns the photospheric driver generating the Alfven 
waves. So far the few observations of Alfven waves in the solar atmosphere seem to indicate 
a broad spectrum range roughly between 100 and 500 s (see introduction) with a power peak 
close to the 5-min power peak of the p-modes. This may indicate that Alfven waves share 
with the p-modes the photosphere as generation site. The waves could be a consequence of 
the constant buffeting of the magnetic field by convection, or possibly, photospheric magnetic 
reconnection. If, instead, Alfven waves are primarily generated by magnetic reconnection 
higher up in the atmosphere it is difficult to explain the observed common power peak. A 
possible scenario coul d be one in whi c h p-m odes trigger magnetic reconnection as in the 



recent simulations by iHeggland et al.l (120091 ). and by some mechanism transfer the same 



periodicity to the waves. In this work we assume mainly that Alfven waves are generated 
in the photosphere and have a white noise spectrum, which would issue from the random 
buffeting motion by convection (or maybe photospheric magnetic reconnection). However, 
we also study the influence of the wave period in the resulting corona. Hence we consider as 
well monochromatic wave drivers generating longitudinal modes with 10, 25, 50, 100, 150, 
200 and 300 s period (as one sinusoidal torsional wave creates two longitudinal waves, the 
periods of the Alfve n waves are twice the periods stated here). This analysis was started in 



Antolin et al.l (120081 ). where it was shown that the most effective waves for coronal heating 
had periods between 100 and 150 s. Here we extend the analysis by studying in detail the 
hydrodynamic response of the loop to the different wave periods. The more realistic random 
photospheric driver is considered for the rest of the parameter space study. 

The standard model of the loop to which we apply the different parameters is the same 
as the one considered in Paper 1. That is, a 100 Mm length loop with an apex-to-base area 
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expansion of 1000, a photospheric magnetic field of 2.3 kG and in which the Alfven waves 
are generated by a random photospheric driver creating a white noise frequency spectrum. 



2.2. MHD equations 

The magnetic flux tube is modeled with the local curvilinear coordinates (s, 0, r) where 
s measures distance along the most external poloidal magnetic field line, is the azimuthal 
angle measured around the rotation axis of the flux tube, and r is the radius of the tube. 
We take the 1.5-D approximation, 

— = 0, — = 0, v r = 0, B r = 0, (1) 

where v r and B r are, respectively, the radial components of the velocity and magnetic field 
in the magnetic flux tube. In the considered approximation, conservation of magnetic flux 
defines the value of the poloidal magnetic field as a function of r alone, B s = B (r /r) 2 , 
where Bo is the value of the magnetic field at the photosphere and tq = 200 km is the initial 
radius of the loop. In the photosphere the value of f3 — 8np/(B 2 + BV) (the ratio of gas 
to magnetic pressures) is unity. We assume an inviscid perfectly conducting fully ionized 
plasma. The effects of thermal conduction and radiative cooling are considered. 

The 1.5-D MHD equations are written as follows. The mass conservation equation: 



the s-component of the momentum equation: 



dv s dv s I dp v\dr 1 B^ d 



the 0— component of the momentum equation: 



the induction equation: 



and the energy equation: 
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where 



V = P — T, e 



1 p 



(7) 
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In the above equations ©-© P, P, and e are, respectively, density, pressure and internal 
energy; v s is the poloidal component of the velocity along the external magnetic field line; 

is the toroidal (azimuthal) component of the velocity; B s and are, respectively, the 
poloidal and toroidal components of the magnetic field; ks is the Boltzmann constant and 
7 is the ratio of specific heats for a monatomic gas, taken to be 5/3. The g s is the effective 
gravity along the external poloidal magnetic field line and is given by 



where g & = 2.74 x 10 4 cm s -2 is the gravity at the base of the loop, z is the length along 
the central axis of the loop, and L is the total length of the loop. The C(t, s) in equation 
(Hj) is a term which simulates the torque motions in the photosphere. It is responsible for 
the generation of Alfven waves. For random photospheric perturbations generating a white 
noise spectrum for the Alfven waves we have 



where Hq denotes the pressure scale height at z — 0. We take Ho = 200 km. The terms 
randl(t) and rand2(t) are non correlated functions that output a number randomly dis- 
tributed between and 1, which changes in time, and / is a parameter that determines the 
strength of the torque. For the case of the photospheric driver generating monochromatic 
Alfven waves we have 



where r is the period of the Alfven waves. We consider waves with periods of 20, 50, 100, 
200, 300, 400 and 600 s. As Alfven waves, through mode conversion, generate longitudinal 
modes having half the period, we obtain longitudinal waves with periods of 10, 25, 50, 100, 
150, 200 and 300 s. Since longitudinal waves are more easily detected during observations 
in this paper we prefer referring to the longitudinal waves that issue from the Alfven waves. 
Hence the periods that are mentioned always refer to periods of longitudinal waves, unless 
stated otherwise. 
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In equation (jHJ), k = 9 x erg s _1 K _1 cm -1 is the Spitzer conductivity corresponding 
to a fully ionized plasma. Here S is a heating term maintaining the initial temperature 
distribution of the loop. The radiative losses R(T) are defined as 



R(T) = n e n p Q{T) 



rr 



Q(T) 



where n = n e + n p is the total particle number density (n e and n p are, respectively, the 
electron and proton number densities, and we assume n e = n p = p/m to satisfy plasma 
neutralit y, with m the proton mass) and Q( T) is the radiative loss function for optically thin 
plasmas (iLandini Sz Monsignori Fossilll99oh which is approximate d with analytical f unctions 
of the form Q(T) = yT 1 . We take the same approximation as in iHori et al.l (119971 ; please 
refer to their Table 1). For temperatures below 4 x 10 4 K we assume that the plasma becomes 
optically thick. In this cas e, the radiative losses R can be approximated by R(p) = 4.9 x 10 9 p 
( (Anderson fc Athayill989l ). In equation the heating term S has a constant non-zero value 
which is non-negligible only when the atmosphere becomes optically thick. Its purpose is 
mainly for maintaining the initial temperature distribution of the loop. 



2.3. Initial conditions and numerical code 

In the present model a sub-photospheric region is considered by adding 2 Mm at each 
footpoint of the loop, in which the radius of the loop is kept constant (hence keeping a 
constant magnetic flux). We take the origin s = z = as the top end of this region. 
The loop is assumed to follow hydrostatic pressure balance in the sub-photospheric region 
and in the photosphere up to a height of AHq = 800 km, where Hq is the pressure scale 
height at z — 0. The inclusion of the sub-photospheric region avoids unrealistic density 
oscillations due to the reflection of waves at the boundaries, thus avoiding any influence 
from the boundary conditions on the coronal dynamics. For the rest of the loop, density 
decreas es as p oc h~ A , where h is the height from the base of the loop. This is based on the 
work bv lShibata et all Jl989al lbl) in which the results of 2D MHD simulations of emerging flux 



by Parker instability exhibit such pressure distribution. The initial temperature all along the 
loop is set at T = 10 4 K. The density at the photosphere (z = 0) is set at p — 2.53 x 10~ 7 g 
cm -3 , and, correspondingly, the photospheric pressure is po = 2.09 x 10 5 dyn cm -2 . As the 
plasma (3 parameter is chosen to be unity in the photosphere, the equipartition magnetic 
field strength B eq = y/Srrpo corresponds to B s = 2.3 kG. For loops having the later as the 
photospheric field strength, the corresponding value of the magnetic field at the top of the 
loop is then B s ^ op = 2.3, 3.8 and 4.6 G respectively for loops expansions of 1000, 600 and 
500. 
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In order to correctly resolve the dynamics in the chromosphere and transition region we 
take a fine spatial resolution in the numerical scheme by setting the grid size As = 0.025H = 
5 km up to a height of ~16000 km. Then, the grid size is allowed to increase as Asj +1 = 
1.03Asj until it reaches a size of 20 km. The size is then kept constant up to the apex of the 
loop. We take rigid wall bo undary conditions at the photosphere. The n umerical schemes 
adop t ed are the CIP scheme ( Yabe fc Aok] 1991) and the M OC-CT scheme ( Evans fc Hawley 



19881 ; IStone fc Normanlll992l ). Refer to lKudoh et al.l (119981 ) for details about the application 



of these schemes. The total time of the simulation is 568 minutes. 



3. Parametric space survey 

In paper 1 it was shown that a corona heated by Alfven waves is characterized by large 
dynamics. High speed plasma flows and large up and down shifts of the transition region 
due to the shocks are common features of this heating mechanism. In order to determine the 
average values of thermodynamic variables in the corona the location of the two transition 
regions are tracked down in time. This is carried out for all the figures displaying coronal 
quantities in the present work. 



3.1. Wave dissipation 



In the present model the linearly polarized Alfven waves have basically only one mecha- 
nism available in order to dissipate their energy: mode conversion. Due to the large density 
stratification of the atmosphere Alfven waves grow in amplitude and the problem easily be- 
comes nonlinear. Alfven waves can then convert a part of their energy to the longitudinal 
slow and fast modes during propagation. This conversion happens mainly while propagating 
from the photosphere to the chromosphere (passing from high-/? to low-/? plasma), although 
important conversion processes also happen in the corona, as discussed below. The result- 
ing longitudinal modes subsequently steep iuto shocks and heat the plasmfl This mode 
conversion mechanism has been shown to be an efficient dissipation mechanism for Alfven 
waves by which the corona can be heated for both open and closed magnetic field re gions 



(IHollweg et al.lll982l : iKudoh & Shibatalll999l : Morivasu et al 



also for the generation of the fast and slow solar wind ( Suzuki fc Inutsuka 20061 ). Another 



2004 



Antolin et al l 120081) and 



dissipation mechanism which may be playing an important role in our model, especially when 



2 It should be noted that being essentially a one dimensional simulation the wave vector is in the same 
direction as the magnetic field, and hence fast modes and Alfven waves are degenerated from each other. 
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considering the monochromatic wave driver in the photosphere, is through resonant damping. 
In this heating scenario closed magnetic flux structures like loops serve as resonant cavities 
for Alfven waves. For certain wave frequencies the loop can resonate thus allowing a high 
transmission of these waves into the corona. This heating mechanism is further discussed in 
section 13.3.11 For the mode conversion mechanism the density fluctuation is crucial for the 
efficient dissipation of the waves. As was shown in Paper 1 and other related works (see for 



instance ISuzuki fc Inutsukal 120061 ) the coronae issuing from Alfven wave heating characterize 
by showing large density fluctuations. Mode conversion takes place wherever the sound and 
Alfven speed have similar values (plasma /3 ~ 1 region), which is mainly the chromosphere, 
but mainly when the nonlinear effects are large. Nonlinear effects can be quantified by the 
ratio of azimuthal to longitudinal magnetic field (or velocities) B^jB s . As shown in Paper 
1, this quantity takes values higher than 0.5 sporadically and ubiquitously in the corona 
due to the large density fluctuation caused by the longitudinal shocks. The large density 
fluctuation throughout the loop also causes the deformation of the Alfven wave shape, which 
causes s mall-scale fluctuati ons of the magnetic pressure and leads to the generation of slow 
mode^l ( Moore et al. 1991 ). Another factor playing an important role in the wave dissipa- 
tion is the wave-to- wave interaction. In this closed magnetic field model we have waves being 
generated at both footpoints, and also reflection of waves at both transition regions. We 
thus have collision of waves and shocks at all times in the corona, which greatly enhances 
nonlinearity and shock strengths. This intermittent mode conversion in the corona leads to 
a bursty intensity profile similar to the usual X-ray inten sity profiles of the sola r corona (as 
those shown by Hinode/XKT for instance). This has led iMoriyasu et al.l (120041 ) to propose 
the idea of nanoflares being the signatures of this nonlinear heating mechanism from Alfven 
waves. It has also been shown from statistics of the resulting heating events that a power law 
distribution of the intensities issues with an index steeper than 2 (cf. Paper 1), a sc enario in 
which the bulk of the heating comes from the small nanoflare-like heating events ([Hudson 
199lh . 



3.2. General characteristics 



Figs. [U [6] and [7| show the average in time of several quantities describing the thermo- 
dynamic state of the corona with respect to the photospheric rms azimuthal velocity field 



3 The mean Alfven speed in the corona obtained from heating by 100 s period waves in our model is 300 
km s _1 . This gives a wavelength on the order of one third of the loop size. Such long wavelengths and the 
large density fluctuations lead to the deformation of the wave shape and the subsequent generation of slow 
modes. 
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created by the torsional motions at the footpoints of the loop for each case of our parameter 
study. In panels a of Figs. [3], [5J [6] and [7] we see the expected general tendency of the mean 
coronal temperature to increase with the amplitude of the photospheric velocity perturba- 
tion. This is expected since higher photospheric transversal velocities mean a higher energy 
flux into the loop (the converse is not always true as we will see below), as shown by panels 
c in Figs. El [5], [6] and U\ which show the variation of the mean rms photospheric magnetic 
energy (B 2 ^/Sn) 1 ^ 2 with the photospheric velocity amplitude. Higher energy fluxes from 
the waves generate stronger longitudinal shocks leading to a higher heating rate as long 
as conditions for mode conversion are met (see sectior j3TTT) . This is shown in the e panels 
corresponding to the same figures, where the volumetric heating rate H is plotted (below 
we explain how we calculate this quantity). We see that conditions leading to long term 
average temperatures above 1 MK are not easily met. Only large expansions (above area 
ratio of 600 between photosphere and corona) together with long lengths (above 100 Mm) 
provide the necessary conditions. Also, the photospheric velocity amplitudes must be large 
enough. In the present case coronae with average temperatures (over the entire simulation 
time) above 1 MK are obtained for {v 2 ^) 1 ^ 2 > 1.6 km s _1 . On shorter timescales, temper- 
atures in the corona can easily reach nanoflare-like temperatures above 1 MK, as shown by 
panels g, which shows the maximum temperatures that can be reached with Alfven wave 
heating. Maximum temperatures show the same dependence with the photospheric velocity 
field as for the mean coronal temperature. Maximum temperatures between 2 and 5 MK are 
found for long period waves (defined here as waves with periods above 100 s). This is further 
discussed in section 13.3.11 As was found in Paper 1 the most efficient wave periods for the 
heating are between 100 and 150 s. The 300 s waves also offer an interesting heating scenario. 
Their high power (as seen in the c panel of Fig. [3]) combines with resonant damping to give 
the highest temperatures in the model. This case is discussed in detail in section 13.3.11 

As the strengths of the shocks increase with the amplitude of the photospheric pertur- 
bations the mass flux into the corona increases due to the higher magnetic pressure from the 
waves, thus increasing the density of the corona (b panels). The Alfven speed consequently 
decreases on the one hand, while the sound speed increases on the other due to the higher 
temperatures, leading to a high increase of the plasma j3 parameter (c panels), which is the 
ratio of gas to magnetic pressures (3 = Snp/ (B 2 + BT). Due to the strong shocks permeating 
the corona, the plasma f3 parameter takes values above 0.1 sporadically and ubiquitously in 
the corona, contributing to mode conversion processes. This behavior is also shared by the 
quantity v^/va, which measures the nonlinear effects and therefore is the main contributor 
for mode conversion. This is crucial for the efficient dissipation of the waves in the corona, 
which therefore gets uniformly heated, as discussed in section [5TTI The average value in time 
of this quantity along the loop is plotted in Fig. [2j Panel a corresponds to our base model: 
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a loop of length 100 Mm, photospheric magnetic field of 2.3 kG and area expansion of 1000 
between the photosphere and the corona, in which a white noise spectrum of Alfven waves 
is generated at the footpoints. Panels b, c and d correspond to variations of the base model: 
with a photospheric magnetic field of 1 kG (panel b), an area expansion of 500 (panel c), 
and a length of 60 Mm (panel d). For each model different values of the rms photospheric 
velocity field are plotted. The solid line corresponds to the value of (v 2 ^) 1 ^ 2 for which the 
highest temperatures are obtained in each model. As a general rule in all our models a high 
value of the ratio v^/va in the corona (> 0.3), or correspondingly in the photosphere (> 
0.060 acquires its minimum value. In the sub-photospheric region the ratio increases again 
due to the low Alfven velocity), is necessary in order to obtain a million degree corona. For 
higher amplitudes of the rms photospheric velocity field (dashed lines) the nonlinear effects 
are reduced (except for the case of the short loop with 60 Mm, which we will discuss in 
more detail in section I3.3.4p and the loop enters a regime in which the corona cannot be 
maintained over long periods of time. This regime is discussed in section 13.2.11 

The uniform heating of the corona can be quantified by comparing the obt ained temper 



ature s with the predicted coronal temperatures based on the RTV scaling law (jRosner et al. 



19781 ). This is demonstrated in Fig. [TJ Panels a, b, c and d display the mean temperature 
in the corona with respect to the predicted temperature from the RTV scaling law, and cor- 
respond, respectively, to the photospheric driver study, the photospheric magnetic field, the 
loop expansion study and the loop length study. It can be noticed that in general the RTV 
scaling law fits remarkably well with the obtained coronal values, mainly for short period 
waves (< 100 s, as discussed in section 13". 3. II) and low photospheric velocity amplitudes. The 
error increase with photospheric velocity amplitude is due to the stronger shocks which in- 
put large amounts of chromospheric material in the corona (creating high spicule-like events) 
and therefore disrupts the local energy balance. Since the RTV scaling law i s found to be 



well s atisfied we estimate from it the volumetric heating rate H. As shown in iRosner et al. 



fjl978[ ). we have: 

H = 9.8 x 10y /6 L" 5/6 . (12) 
The estimations resulting from this equation are plotted in the e panels of Figs. [3j [5j [6] and 

El 

Despite the higher momentum of the waves as the amplitudes of the photospheric driver 
increase the upward shifts of the transiti on region (spicule-like events) remain roughly con- 
stant. This has been previously found by lHollweg et al. Jl982h . As explained in that paper 



for the case of open flux tubes, given a high amplitude shock we will have a faster upward rise 



4 The photosphere is pointed out in the figure by the location close to the footpoints where |u^|/va 
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of the transition region and thus also a longer duration for the following shock to reach the 
transition region, hence giving the latter more time to recede from gravity to lower heights. 
In the present case of closed flux tubes we have found the same behavior, thus producing a 
mean length of the corona that is roughly constant as seen in the f panels. However, when 
the photospheric velocity amplitude exceeds a critical value this is no longer the case. In- 
deed, a sharp decrease of the coronal length can be noticed when the photospheric azimuthal 
velocities exceed a certain critical value which depends mostly on the geometrical parameters 
of the loop (length and expansion). The critical photospheric velocity decreases the longer 
the periods of the waves as discussed in section [3.3. II For photospheric velocities higher than 
the critical velocity the corona enters an unstable regime in which the temperature, density 
(and consequently the plasma j3 parameter) stop increasing. These cases are marked with 
blank symbols in all panels. We will now discuss this corona collapse regime. 

3.2.1. Corona collapse 

In order to create a corona wave heating coming from mode conversion of Alfven waves 
(cf. section I3.ip has to exceed the losses due to radiation and (downward) conduction. 
The rapidity with which this energy balance is reached and a corona forms will depend on 
many parameters, such as the initial conditions as well as geometrical aspects (length and 
expansion of the loop for instance). In panels h of Figs. [31 El M arid [7] we show the time span 
of the corona over the entire simulation time (568 min). We can see that for each set of the 
considered loop parameters there is a range of photospheric velocity amplitudes for which 
the energy balance in the corona is maintained over most of the simulation, i.e. for which 
the corona is stable. The time span of the corona in this region is roughly constant. However 
there is a critical velocity above which the stability of the corona is compromised and the 
time span drops off sharply. These cases are marked with blank symbols in all panels. For 
these cases the corona can break down at some point and reform at a later moment, or not be 
reformed at all. More specifically, we have defined a stable corona as a corona which lasts for 
more than 60 % of the total simulation time, which corresponds roughly to 5.7 hours. The 
value of the photospheric magnetic field (and hence the value of the background magnetic 
field) does not seem to influence the critical photospheric velocity at which the instability 
sets in, as seen in the h panel of Fig. [51 For all the values considered for the photospheric 
magnetic field the corona gets unstable for (v 2 ^) 1 ^ 2 > 2.2 km s _1 (see section 13". 3. 2p . Loops 
with thinner area cross-sections and shorter lengths have much smaller stable regimes (see 
sections 13.3.31 and 13.3.41) . Unstable coronae are seen to happen only for wave periods above 
100 s. The critical velocity amplitude above which instability sets in is also dependent on 
wave period. The longer the period the lower the critical velocity. For instance, for 100 and 
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200 s waves we find unstable coronae already at 2.8 km s 1 and 1.8 km s 1 respectively. 

The unstable regime can be understood partly by analyzing the behavior of the volu- 
metric heating rate. As shown by the e panels of Figs. El El El and [7] this quantity reaches 
a maximum for values of the photospheric velocity field which are generally lower than the 
critical velocity for which the corona collapses. The volumetric heating rate is then either 
constant or slightly decreases. Correspondingly, the mean temperature and density in the 
loop reach maximum values after which they are either constant or decrease. The presence 
of these maxima can be understood by means of the energy balance in the loop, which can 
be approximated in the following way: 



n 



2 



T 7/2 



H=-Q(T) + k — , (13) 

where the first term on the right corresponds to the radiative losses, and the second term cor- 
responds to thermal conduction. For values of the photospheric velocity field lower than the 
critical value we have a constant volumetric heating rate H max , temperatures corresponding 
to an optically thin cooling which can be approximated by Q(T) = yT" 1 / 2 with y = 1Q ~ 18 - 48 



(we take the same approximation for temperatures above 10 51 as in lRosner et al.lll978l ) and 
a roughly constant length of the corona L. The maximum value of the volumetric heating 
rate then naturally sets a maximum value for the density and temperature. We have: 

4 , T 7 ' 2 , 

n = [-T l '\H max - no—)} 1 ' 2 (14) 
X £ 

This equation defines a maximum density value n max = \^H max TU' 1 ^ L ^ 2 for a temperature 

T m = { Hv l a *^ 2 ) 2/7 ■ Replacing by numerical values: H max ~ 3 x 10~ 5 erg cm" 3 s -1 , L ~ 8 x 10 9 
cm, and kq ~ 10~ 6 erg s" 1 K" 1 cm" 1 , we obtain a maximum density of n max ~5x 10 8 cm -3 
for a temperature of T m ~ 6 x 10 5 K, which roughly match the critical temperature and 
density values we obtain. 

Fig. [2] shows that nonlinear effects are reduced in the unstable regime. Hence, mode 
conversion is reduced leading to a saturation of the volumetric heating rate. The saturation 
is also due to the height of the spicules. As explained in the previous section, the height of 
the transition region is roughly unaffected by the amplitude of the photospheric driver, as 
long as this one does not exceed the critical value. When the critical amplitude is exceeded 
spicules become much longer due to the decreasing effective gravity along the loop. The 
cooling effect on the corona is therefore large and the volumetric heating rate from the 
waves is saturated. Correspondingly, for values of the photospheric velocity close to the 
critical velocity the energy balance along the loop is not well maintained, and the error 
from the RTV scaling law increases, as seen in Fig. [IJ This results in the collapse of the 
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corona. As stated previously, long period waves have shorter critical velocity than short 
period waves. This is due to the stronger and less frequent shocks obtained with longer 
periods as discussed in section 13.3.11 For 300 s waves the obtained coronae are unstable even 
with low amplitude photospheric velocities (except for only one case). Long period waves 
can carry large amounts of chromospheric material into the corona, increasing their density 
as shown by panel b of Fig. [31 Since the energy flux into the corona is constant in time the 
heating per unit mass decreases in time, an effect which is stronger for long periods due to 
the lower number of shocks in the corona at a given time. For low amplitudes of the long 
period waves (that is, as long as the cooling effect from the spicules is not important) the 
coronae reach a critical state in which they are susceptible to a thermal instability, by which 
the corona rapidly cools down (in the order of minutes) and the loop is evac uated. A similar 



cooling event is obtained for loops which are heated tow ards their footpoints (lAntiochos et al. 



19991 ; iMiiller et al.ll2003l ; iMendoza-Briceno et al.ll2005l ). a phenomenon termed "catastrophic 



cooling", which has been proposed as the physical explanation for coro nal rain. The cyclic 
pattern of heating and cooling in the loops ("limit cycles", as termed by lMiiller et al.ll2003l ) 
is only seen for low amplitudes of the monochromatic driver generating long period waves. 
In all other cases, for instance for thinner and shorter loops or for high amplitudes of the 
random photospheric driver, the collapse of the corona that is obtained is of a different nature 
than that of catastrophic cooling. In these cases the collapse of the corona is not due to a 
local increase of the radiative losses in the corona from the formation of a cool condensation 
as in the case of catastrophic cooling. Instead, the long spicules shorten dramatically the 
length of the corona, disrupting the local energy balance, and no condensation is formed. 
Whether Alfven waves can reproduce coronal rain, as catastrophic cooling does, despite the 
uniform heating they produce is an interesting question which will be further investigated 
in a following paper. We will now discuss individual features of each set of parameters. 



3.3. Individual characteristics 

3.3.1. Photospheric driver 

The loop model in this section consists on a 100 Mm length loop with an expansion 
factor of 1000 and a photospheric magnetic field of 2.3 kG, and in which monochromatic 
Alfven waves are generated at both footpoints. 

Panels a to h in Fig. [3] display the characteristics of the coronae produced by monochro- 
matic slow and fast MHD waves resulting from the mode conversion of Alfven waves. We 
can see that short period (^ 50 s) and long period (^ 100 s) waves produce very different 
coronae. Long period waves produce hotter and denser coronae than short period waves as 
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shown by panels a and b. Coronae produced by long period waves are subject however to 
the unstable regime discussed in section [3.2.11 when the photospheric velocity amplitudes are 
high, as can be seen in panel h, and by the blank symbols in all the panels. For wave periods 
below 150 s, for the same velocity amplitude we obtain higher coronal temperatures the 
longer the period, and we find temperatures above 1 MK for velocities (v 2 ) 1 ^ 2 in the range 
0.9 — 1.6 km s" 1 for 100 and 150 s period waves. Waves with lower periods (high frequency) 
are rapidly damped when propagating through the lower atmospheric layers. On the other 
hand, since longer periods imply longer wavelengths the number of shocks heating the corona 
at a given time is reduced and the corona is less uniformly heate d. This is shown in panel a 
of Fig. [TJ, where the match with the RTV scaling law is plotted (IRosner et al.lll978l ). Long 
period waves see m to suffer also from a period cut-off, reducing their efficiency for heating 
the solar corona. iMusielak et ajj (120071 ) have shown that torsional Alfven waves have cut-off 
free propagation for the c ase of thin magnetic flux tubes, an d thus can carry large amounts 



of energy into the coro na. IVerth et al. 



2009. in preparation! ) have pointed out however that 



the assertion made by IMusielak et al.1 (120071 ) is valid only when the temperature in the flux 



tube does not differ from that of the external plasma. In the present case the flux tube 
presents velocity gradients (which issue from the densit y stratification and expansion of the 
magnetic flux tube) which introduce a period cut-off. iRouth et al.1 (120071 ) have considered 
this problem analytically and shown that the cut-off period for torsional Alfven waves is local 
(as opposed to global, where coupling to gravity or density stratification exists), introduced 
by gradients of the characteristic wave velocities. Caution should be taken however when 
applying the results of the previous paper, sinc e it contains an approximation whi ch may not 
be valid for thick flux tubes, as pointed out by lVerth et al.1 (12009. in preparation! ) . However, 
their main result concerning the existence of the cut-off period should remain valid. Hence, 
we can expect the cut-off period to vary with height in the same way as the velocity of the 
torsional waves va varies. In view of the questionable approximation taken in their paper, 
we can take their reported cut-off period as an approximation. In the case of propagation 
through an isothermal medium they have shown that the period cut-off is given by 



PcuM = (15) 
va{s) 

where they approximate va(s) = u^ i0 exp(s/m/i), with s the propagating distance, h the 
pressure scale height and m is a positive scaling factor. According to this result the propa- 
gation of Alfven waves into the corona in our model should be significantly reduced. In our 
model the exponential increase in the Alfven velocity happens in the upper chromosphere. 
In that region a typical value for the mean sound speed is c s = 15 km s _1 , giving a pressure 
scale height of h ~ 500 km, and the mean Alven velocity is va = 10 km s -1 . Hence we get a 
cut-off period of ~ 50 s for the torsional Alfven waves, so 25 s for the longitudinal waves. This 
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implies that effects from the cut-off period should be present in most of the cases considered. 
Since propagation of the Alfven waves do happen (at least to heights where mode conversion 
occurs) resulting in the formation of a corona, the discrepancy could be due to the highly 
nonlinear effects in this model (and also to the questionable approximation of the paper). 
The importance of the nonlinear effects for Alfven wave heating is thus emphasized. Apart 
from the cut-off period for Alfven waves the normal cut-off period for longitudinal waves is 
also present. In our model, the value of this period cut-off is close to 200 s. Hence, a large 
fraction of longitudinal waves with periods larger than 200 s issuing from mode conversion 
occurring below the transition region won't be able to propagate into the corona. This may 
explain the low heating efficiency from the 200 s period waves. For the 300 s waves however, 
we still can have a hot corona as seen in the a panel. We will discuss why this is so further 
ahead. 

Long period waves can create very hot heating events close to 5 MK as show n in panel g, 



accou nting for the observed nanoflares with instruments such as Hinode/XRT (IGolub et al. 



20071 ). Despite the fact that short period waves do not contribute substantially to the hot 
temperatures of the corona they do have an important effect for the thermal stability of the 
loop. As previously stated the coronae obtained with short period waves are more uniform 
due to the higher number of shocks permeating the corona. The large cool mass upflows into 
the corona due to the waves are reduced, resulting in longer lengths and time spans of the 
corona (panels f and h). 

In order to better understand these results, we plot in panel c of Fig. [3] the mean az- 
imuthal magnetic energy B^/8tc at the photosphere with respect to the photospheric velocity 
amplitude for each period. is a measure of the twist exerted at the footpoints of the loop, 
which increases with the period of the wave. Indeed, we can see that for a same input of 
magnetic energy, there is a shift towards lower amplitudes of the photospheric velocity field 
the longer the period of the wave. Hence, it is to be expected that longer period waves 
produce hotter coronae as shown by panel a. Also, larger twists imply higher wave pressure 
and mass flux into the corona, thus explaining the higher densities obtained with long period 
waves. The 300 s waves exhibit a considerable shift: with a photospheric velocity amplitude 
as low as 0.25 km s _1 , the corresponding magnetic energy input is as high as for the 100 s 
period waves with an amplitude of (f^) 1 ^ 2 = 1.3 km s _1 (namely, ~ 1.8 x 10 4 erg cm). This 
explains why we obtain the same coronal temperatures for both cases (panel a). 

The large difference that can be seen in the photospheric azimuthal magnetic energy 
(panel c), and in general in all panels of Fig. [3] between the 300 s waves case and the other 
periods points to the presence of an additional heating mechanism. Indeed, apart from mode 
conversion there is another heating mechanism which seems to be playing an important role, 
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especially when considering a monochromatic wave driver. Closed magnetic structures such 
as loops may act as resonant cavities for the Alfven waves, which can suffer from reflection 
due to the steep gradients of the Alfven velocity and be trapped in the corona. Resonances 
can increase the transmission coefficient from the waves, thus allow ing large ener gy fluxes 



into the corona. Heating through resonances was first considered by llonsonl (119781 ). 



In Fig. H] we plot the map of the power spectrum at each point along the first 20 
Mm of a loop with 300 s monochromatic waves for the longitudinal velocity component v s 
corresponding to the longitudinal modes, and azimuthal velocity component corresponding 
to the Alfven waves. In this case the rms photospheric velocity amplitude is < > 1//2 ~ 0.75 
km s _1 , and we can see a clear resonance pattern from the corona down to the chromosphere. 
In the longitudinal velocity map, the fundamental mode at 300 s and the first (150 s), second 
(100 s) and third (75 s) harmonics coming from the mode conversion of the Alfven waves 
can be clearly seen. However, in the azimuthal velocity map the fundamental mode at 600 
s and only the even harmonics (second, fourth, and so on) can be seen. Since we have 
set the monochromatic driver to generate waves that are in phase, the absence of the odd 
harmonics is simply due to a destructive interference between the Alfven waves generated at 
both footpoints. Nevertheless, since mode conversion can happen in a time scale of minutes, 
before the disappearance of the odd harmonics mode conversion from these waves occurs 
resulting in the odd harmonics of the longitudinal velocity map. 



Hollwegi (11984 ) was the first to consider resonant damping in the case of a loop, where 
the later is approximated by a three-layer model in which the Alfven velocity increases 
exponentially i n the bo u ndary layers (chromosphere and transition region) and is constant 



in the corona. iHollwegi (11984J ) showed that when a wave is generated at one footpoint the 



transmission coefficient into the corona can be written 

rM|1+( v^>) , _ (16) 

where va, c is the Alfven speed in the corona, L c is the length of the corona, h is the pressure 
scale height, u is the frequency of the wave and k = u/v^c is the wavenumber. The coronal 
resonances are excited for the periods 

(17) 



mv Ac ' 



where m = 1,2,3.... At these periods the transmission coefficient has its maximum value. 
Eq. [171 is simply the travel time for the Alfven waves back and forth between the boundaries of 
the resonant cavity, and hence denotes the normal oscillating modes for the Alfven waves. In 
that work the height of the transition region was however not considered, which is important 
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in our model. Indeed, we have seen that the height of the spicules can be large. Furthermore, 
we also consider a sub-photospheric region at both footpoints. We can estimate the resonant 
periods in our model by adding the travel time of the waves in the regions below the transition 
regions. A more precise estimation of the resonant periods would take us too far from our 
present parameter study, and thus is left as future work. We have: 

2L ?m 2L r 

r res = 2 ^ + c -, (18) 

mvAjp mv A ,c 

where Lf p is the length of the region below the transition region, and v aj p is the mean Alfven 
velocity in that region (all quantities denote averages over time). Inserting mean values of 
the 300 s waves case where resonance is detected: L c = 84 Mm, Lf p — 10 Mm, va, c = 440 
km s _1 and v A ,f P = 190 km s -1 , we obtain: P ~ 590/m s, which matches well the resonant 
periods of the Alfven waves (and hence the resulting periods for the longitudinal waves). 



3.3.2. Photospheric magnetic field 

In Fig. [5] we plot the results of the parameter survey for different values of the photo- 
spheric magnetic field B S>Q . We test B s Q = 2.3 kG, 1.5 kG and 1 kG. Here the model consists 
of a 100 Mm length loop with an expansion factor of 1000 and a random photospheric driver. 
Having set in our model the value of the photospheric pressure to po = 2.09 x 10 5 dyn cm -2 , 
the equipartition magnetic field strength B eq = \^8np^ corresponds to B s0 = 2.3 kG. In 
this case we use the more realistic random photospheric driver for the generation of Alfven 
waves. From the panels of Fig. [5] we can see that the thermodynamical properties of the 
corona do not change considerably with the value of the photospheric magnetic field. For 
low amplitudes of the photospheric velocity field ((v 2 ) 1 ^ 2 < 1.1 km s _1 ) the photospheric 
magnetic energies, temperatures, densities, lengths and time spans of the corona are similar 
as shown by panels a, b, c, f, g and h. Also, the critical velocity above which the corona 
is unstable is roughly the same: (w 2 ,) 1 / 2 ~ 2.2 km s" 1 . The a and b panels in Fig. [2] show 
that the nonlinear effects are also roughly similar, since these are rather independent on the 
background magnetic field. For photospheric velocities between 1.2 and 2.2 km s _1 , however, 
hotter and denser coronae are obtained for stronger photospheric magnetic field, as shown 
by panels a and b in Fig. [5J For instance, for a photospheric velocity amplitude of ~ 2 km 
s _1 the coronal temperature and density are, respectively, 7.4 x 10 5 K and 2.5 x 10 8 cm -3 
for B St0 =1 kG, and 1.1 x 10 6 K and 5.5 x 10 8 cm -3 for B St0 = 2.3 kG. However, the plasma 
j3 parameter is roughly equal for velocities (v^) 1 / 2 > 2 km s _1 as seen in panel d. For lower 
values of the photospheric velocity field we obtain higher plasma (3 the lower the value of 
the photospheric magnetic field. 
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The disparities between the three cases can be understood by the different pressure 
distributions along the loop obtained by varying the photospheric magnetic field. Since 
the equipartition magnetic field corresponds to the case B s0 = 2.3 kG, the other cases 
have higher plasma (3 at the same height in the photosphere and the chromosphere. As 
the amplitude of the torsional motions in the photosphere increases the magnetic pressure 
increases and pushes photospheric and chromospheric matter upwards. Since the magnetic 
pressure depends not only on the square of the torsional motions but also on the square of 
the longitudinal magnetic field the uplift effect is stronger for stronger photospheric magnetic 
field. Then, the gas pressure p increases faster the higher the magnetic field. A higher coronal 
gas pressure implies a higher volumetric heating rate as shown by Eq. [12], and illustrated by 
the e panel of Fig. [51 thus explaining the higher temperatures. 

New observations with the Solar Optical Telescope on board of the Hinode satellite have 
shown patches of strong magnetic of 1 - 1.5 kG sc attered uniformly all over the surface of the 



Sun, especially in Quiet Sun and Polar regions (ITsuneta et al.ll2008l ). Our results indicate 
that Alfven waves that may be generated at the footpoints of loops spanning from these 
regions do not have energy limitations as far as the photospheric magnetic field is concerned. 



3.3.3. Loop expansion 

We now consider different expansions of the loop displaying area ratios of 1000, 600 
and 500 between the corona and the photosphere. The loop is 100 Mm in length and has 
a photospheric magnetic field of 2.3 kG. The Alfven waves are generated with the random 
photospheric driver. 

Fig. [6] displays the influence of the expansion of the loop in the thermodynamic proper- 
ties of the corona. As the expansion of the loop is reduced the temperature and density in 
the corona are reduced as well, as shown by panels a and b, despite the same photospheric 
magnetic energy at the footpoints (panel c). For the loop with an expansion factor of 500 
all the average temperatures in the corona are well below 1 MK (in this case the maximum 
achieved average temperature for a stable corona is 7.4 x 10 5 K, which appears to be an 
exception due to the large difference with the other cases). Consequently the plasma (3 
differs roughly by one order of magnitude in the corona between the loops with expansion 
factor 1000 and 500 (panel d). Together with the temperatures and densities of the corona, 
the critical velocities above which we have unstable coronae are reduced as well with lower 
loop expansions, as shown by panel h. While the critical velocity above which the corona 
collapses is (v^) 1 ^ 2 ~ 2.2 km s _1 for a loop expansion of 1000, the corresponding critical 
velocities for loop expansions of 600 and 500 are ~ 1.7 and ~ 1.3 km s _1 respectively. In the 
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stable velocity range the lengths of the coronae are comparable, although there are 2 cases 
for the expansion factor of 500 for which the corona is barely stable and the loop lengths are 
reduced by half, namely, (f|) 1//2 ~ 0.9 and ~ 1 km s _1 . We have also performed simulations 
with loop expansions of 300 but no stable corona was found. 



As discussed in section 12. 1[ the expansion of magnetic flux tubes is a parameter which 
can greatly constrain the wave heating. Thin flux tubes allow large energy flows from 
Alfven waves into the corona since they do not introduce any Alfven velocity gradients (as 
long as the external temperature and internal temperatu re of the flux tubes do not differ), 



thus allowing cut-off free propagation of Alfven waves (IMusielak et al.l 120071 ; IVerth et al. 



2009. in preparation! ) . When expansion is strong enough to introduc e velocity gradient s the 
spectrum of Alfven waves that do not become evanescent is reduced (IRouth et al.ll2007l ) , as 
discussed in section 13.3.11 (caution s hould be taken however when app lying the results stated 
in the later paper, as pointed out by lVerth et al.ll2009. in preparation! ) . However, in our case 
we have lack of heating for thinner flux tubes, while the strong expansion case of 1000 can 
easily achieve temperatures above 1 MK. 

A simple calculation can illustrate the dependence of mode conversion on the expansion 
of the loop. Assuming a typical Alven speed of 10 km s _1 in the upper photosphere/ lower 
chromosphere and a period of 100 s we obtain a wavelength of 1000 km for Alfven waves 
propagating from the photosphere. Now, all the loops considered here display roughly the 
same initial expansion in the photosphere. We have an expansion of roughly 1000 km in the 
first 1000 km above the photosphere, thus m atching the wave length of the Alfven wave. A 
slightly higher expansion has been reported by I Jess et al.l (120091 ) by comparing Ha continuum 
with Ha core images in observations of a bright point with the SOUP instrument of SST. 
Since this region is the region where the plasma (3 is close to unity we have ideal conditions 
for mode conversion for the Alfven waves. Fo r instance, deform ation of the wave shape 
and excitation of slow modes should take place (IMoore et al.lll99ll ). Since our 3 cases have 
roughly the same initial expansion, this implies that mode conversion taking place higher up 
in the chromosphere is very important for the Alfven wave heating to be effective. 

The most important influence on the heating from the expansion of loops in our case 
has to do with the nonlinear effects. Higher expansions produce higher azimuthal velocities 
in the corona, which influence flows along the loop through the centrifugal force (cf. third 
term on the right hand side of equation [3]). This in turn increases the density in the corona 
as seen in panel b of Fig. [6j thus decreasing the Alfven velocity. Consequently the ratio 
v <j>l v Ai which measures the nonlinear effects, increases. This is illustrated in panels a and 
c of Fig. [21 where this ratio is plotted for loop expansions of 1000 and 500, for different 
amplitudes of the photospheric velocity field. From basic concepts it can be easily shown 
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that 

^ oc V /4 , (19) 

VA 

where A is the area expansion coefficient. As seen in Paper 1, v^/v s becomes high ubiqui- 
tously and sporadically in the corona. When this happens nonlinear effects are large and 
energy is transfered from the Alfven mode to the longitudinal slow and fast modes, which 
steepen into shocks and heat the plasma (see section IBTTj) . Hence, the reason for the absence 
of heating in our thin loops is not the lack of energy flux or larger radiative losses but a lack 
of dissipation of the waves. 



3.3.4- Loop length 

Observations of coronal loops and posterior analytical studies based on the energy bal- 
ance in these closed magnetic field structures have shown that the coronal temperature in uni- 



form ly heated loops increases with the length of the corona (IRosner et al.lll978l ; ISerio et al. 



198ll ). In order to study the effect of length on Alfven wave heating here we consider coronal 
loops with lengths of 100, 80 and 60 Mm. The loops have the same expansion factor, namely 
1000, and a photospheric magnetic field of 2.3 kG. The waves are generated with the random 
photospheric driver. 

Fig. [7] shows the effect of loop length in the creation of the corona from Alfven wave 
heating. The results in this case are similar to the previous case in which the expansion of 
the loop is considered. We can see that as the length of the loop is reduced the heating of 
the loop is dramatically reduced (panel a), despite the similar magnetic energy input at the 
footpoints (panel c). Already for the 80 Mm length loop we cannot obtain a corona with an 
average temperature above 1 MK. Densities are also reduced and therefore also the plasma (3 
(panels b and d). The sensitivity of the coronae to the loop length can also be seen in panels 
f and h, where we can see that the range of photospheric velocities over which the coronae 
can be created and maintained is considerably reduced with respect to all the previous cases. 
The critical velocities above which the coronae collapse are (u? y/ 2 ~ 1.1 and ~ 0.6 km s _1 , 
respectively for 80 and 60 Mm lengths. For photospheric velocity ranges in which the corona 
is stable we can see from panel f that the lengths of the corona vary considerably. For L = 80 
Mm the lengths vary between 25 and 50 Mm. For L = 60 Mm we have only 2 stable cases 
with lengths between 10 and 20 Mm. 

The lack of wave heating can be understood partly by considering the number of shocks 
at any time heating the corona. In section 13.3.11 we showed that long period waves produce 
hot coronae while short period waves account for their stability by heating them uniformly. 
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The most efficient period range for the heating was found to be between 100 and 150 s for 
the longitudinal waves, thus 200 and 300 s for the Alfven waves. For a photospheric velocity 
of ~ 1 km s _1 typical Alfven velocities in the corona are ~ 400 and ~ 300 km s _1 for loop 
lengths of 80 and 60 Mm respectively (as derived from panel b). Taking a period of 200 
s, this gives wavelengths of 80 Mm and 60 Mm for the Alfven waves respectively, which is 
just the length of the loops. We see that the wave has barely the distance to propagate one 
wavelength before reaching the other footpoint, where the dense sub-photospheric regions 
will not allow an easy reflection. Taking a roughly constant sound speed in the corona it 
is easy to see that the necessary distance for shock formation will be at least on the order 
of a wavelength (in the case of a high velocity amplitude). Since both footpoints generate 
Alfven waves we then may have a maximum of 4 shocks propagating in the corona at any 
time (neglecting all other dissipative mechanisms such as wave-to- wave interaction), and 
may not easily reflect in the transition regions due to their long wavelengths. The same 
reasoning applied to a 100 Mm length loop gives roughly 12 shocks propagating at any 
time. Furthermore, in the later case shocks can easily reflect at the transition region and get 
trapped in the corona, thus increasing the dissipation in the corona. 

Fig. [2] shows that nonlinear effects are considerably reduced for shorter loop lengths. 
The highest value of the ratio v^/va for stable coronae is about 0.1, which is not enough for 
mode conversion to be efficient. The main reason for the lack of hot coronae in short loops 
is the same reason for their instability. Due to the shorter length the effective gravity is 
significantly reduced and the height of the spicules increases considerably as shown in panel 
f of Fig. [7J As the amplitudes of the waves increase large amounts of cool matter are input 
into the loop. Spicules have then a stronger cooling effect the shorter the loop length. 



Discussion and conclusions 



Alfven waves constitute one of the main candidates for heating the solar corona. The- 
oretically easy to be generated from convective turbulent motions in the photosphere it has 
been shown that these waves are able to carry sufficient amounts of energy in order to bal- 
ance energy loss es from radiation and conduction and hea t the corona up to the observed few 



million degrees (Uchida fc Kaburakil 



1989; Ruderman et al. 



1997 



1974 



WentzejJl97J ; iHollweg et al.lll982l : IPoedts et al. 



Kudoh fc Shibatalll999l ). Finding a suitable dissipation mech- 



anism is a very difficult task, which has spawned a large active research community in solar 
physics. In this work we have considered a model in which mode conversion acts as the dis- 
sipative mechanism. More precisely, through nonlinear effects based on density fluctuations 
and wave-to-wave interaction in the chromosphere and corona the Alfven mode is able to 
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transfer some part of its energy to the fast and slow modes, which steepen into shocks and 
heat the plasma. 

We have conducted a parameter survey in which the effect of geometrical quantities 
such as the loop expansion and the loop's length on the efficiency of Alfven wave heating is 
studied. We further investigate the effect on the heating of other physical parameters such 
as the photospheric magnetic field and the generation of Alfven waves with a monochromatic 
and white noise spectrum in the photosphere. 

In Fig. [8] we plot the volumetric heating rate with respect to the considered parameters 
of the loop: period of the monochromatic longitudinal waves issuing from mode conversion 
of Alfven waves (top left panel), photospheric magnetic field (top right panel), expansion of 
the loop (bottom left panel), and loop length (bottom right panel). The volumetric heating 
rate for each case corresponds to its extremum values in the range between (v^) 1//2 ~ 1.3 and 
2.3 km s -1 found in the e panels of Figs. [3j [5], [6] and which is the most efficient amplitude 
range for coronal heating we have found. We can see that the volumetric heating rate is an 
increasing function for all the considered parameters in the chosen azimuthal photospheric 
velocity interval. This figure summarizes well the results of our parameter survey. As the 
amplitudes of the twists in the photosphere generating the Alfven waves increase the momen- 
tum and energy flux of the waves increase. We have found that nonlinear effects generally 
increase as well, and mode conversion happens not only in the chromosphere but ubiqui- 
tously in the corona, thus heating the plasma. Suitable conditions for Alfven wave heating 
are thus strongly dependent on the nonlinear effects. The importance over nonlinearity fa- 
vors thick and long loops (bottom left and right panels), and strong photospheric magnetic 
fields (top right panel). The dependence on the later is however not crucial and flux tubes 
with 1 kG field concentrations at their footpoints may be heated efficiently by Alfven waves. 
On the other hand, the expansion and the length of the loop are crucial parameters for the 
heating. As we can see in Fig. [HI the volumetric heating rate in the bottom left and right 
panels does not vary much respect to the loop expansion and the length, which shows the 
high sensitivity of the thermodynamics in the loop on the heating. Nonlinearity is directly 
proportional to the loop expansion. Furthermore, the strong expansion of loops from the 
photosphere to the transition region happens in a height comparable to the wavelength of the 
Alfven wave, which leads to the deformation of the Alfven wave and excitation of slow modes 



( iMoore et al.lll99ll ). Since the plasma (3 parameter in that region is close to unity, further 
mode conversion is expected. Alfven waves have long wavelengths and in order for them to 
convert into longitudinal waves which steep into shocks they need to propagate relatively 
long distances. Loops with lengths lower than 80 Mm cannot be heated by Alfven waves in 
the present model. Caution must be taken however when drawing conclusions from these 
results. When calculating the lengths of loops in observations what is actually calculated 
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is the length of the corona, which is the visible part of loops in EUV or X-ray images. In 
our model, the lengths of the coronae with average temperatures of ~ 1 MK we obtain with 
Alfven wave heating oscillate between 50 and 80 Mm with a mean around 70 Mm. Further- 
more, we have defined a stable corona, as a corona which can be maintained over roughly 
5.7 hours. However, observations show that coronal loops are dynamical entities which ex- 
hibit heating and cooling processes constantly. Likewise, the loops with a collapsing corona 
obtained here exhibit heating events with nanoflare-like temperatures constantly, which can 
match the observed bursty X-ray intensity profiles and the statistics displaying power law 
distribution of the heating events (cf. Paper 1). Also, a loop which is visible in soft X-rays 
over a large period of time may in fact result from the effect of different threads in the 
loop being heated at differe nt times, thus giving the impres sion of a steady and uniformly 
heated loop (as discussed in iPatsourakos fc Klimchukll2009l ). In the simple model we have 
considered we have assumed that the heating acts uniformly along the radial direction, thus 
disabling the possibility of many threads being heated at different times. We thus cannot 
completely reject the possibility of Alfven waves heating the coronae of short and thin loops. 
This idea should however be tested in an extended (at least 2.5 dimensional) version of this 
model. 

All coronae show a range of velocities for which they are stable throughout the simu- 
lation and for which their lengths are constant. However, when the photospheric velocity 
field exceeds a critical value the corona collapses. The critical velocity is dependent on the 
parameters of the model. In this regime wave heating can no longer account for the large cool 
mass flux from the increasing wave amplitudes and collapses. Generally, the collapse is not 
due to a local increase in the corona of the radi ative losses from the fo rmation of a cool con- 
densation as in the case of catastr ophic cooling ( Antiochos et al.|[l999l ). and does not exhibit 
"limit cycles" (jMiiller et al.l 120031 ). The later are characteristic of the catastrophic cooling 
mechanism which is proposed as an explanation for coronal rain. Further investigation in 
that direction is however required and will be the subject of a future paper. 

All of the obtained coronae are uniformly heated irrespective of the parameters in the 
model. This is a characteristic of Alfven wave heating. The distribution of the heating in 
coronal loops as derived from observations is a matter of debate. Indeed, different methods 
for the analysis of observational data may lead to different results. A famous example of 
this fact is the analysis of the data set from soft X-ray observations with the Yohkoh/SXT 
instrument. An interpretat ion of the h eating distribution of the observ ed loops in terms o f 
apex concentrated hea ting jRealell2002h. footpoint concentrated heating (jAschwandenll2001bl ) 
and uniform heating (IPriest et al.lll998l ) has been given. There seems to be however more 
observational evidence of coronal loops being heated towards the footpoints in active regions 
( lAschwandenl l2001al ; lAschwanden et al.l 120001 ) . Further evidence of this fact has been found 
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by lHara et al.l (120081 ) using the Hinode/ELS instrument, in which active region loops are 
shown to exhibit upflow motions and enhanced nonthermal velocities in the hot lines of Fe 
XIV 274 and Fe XV 284. Possible unresolved high-speed upflows were also found. In Paper 
1 we found that footpoint or uniform nanoflare heating exhibit hot upflows, thus fitting in 
the observational scenario of active regions, while Alfven wave heating was foun d to exhibit 



hot downflows, which may fit in the observational scenario of quiet Sun regions (IChae et al. 



19981 ; iBrosius et al.l 120071 ). The uniform heating from Alfven waves further supports this 
conclusion, since it is a characteristic of loops more generally found in quiet Sun regions. 
Active region loops exhibit low expansion factors due to the high magnetic field filling factor 
in those regions. In this chapter we have found that Alfven wave heating is effective only in 
thick loops (with area expansions between photosphere and corona higher than 600), further 
emphasizing our conclusion that active regions may not be heated by Alfven waves. These 
results seem to point to an important role of Alfven wave heating in Quiet Sun regions, 
where loops are often long, expand more than in active regions, and kG (or higher) bright 
points are ubiquitous. 

Another important result we have found is that both long (above 100 s) and short period 
(below 50 s) waves (where the period is the resulting period of the longitudinal modes from 
mode conversion; the period of the Alfven waves is twice the stated period) play an important 
role in the heating of the loops. Long period waves produce very hot heating events in the 
corona and thus increase the average temperature of the corona due to the strong shocks 
they produce. Short period waves are responsible for keeping the corona, thus making it 
stable through uniform heating from the numerous weak shocks they produce. We have 
therefore a compromise between long and short period waves leading to efficient Alfven wave 
heating. The 300 s period waves are found to have considerable power to heat the corona 
even with low photospheric velocity amplitudes of ~ 0.2 km s _1 . For these waves resonances 
seem to be triggered in the loop, which allow higher transmission from the waves into the 
corona and subsequent mode trapping leading to efficient dissipation. This resonant damping 
mechanism is however not observed for random wave generation (white noise spectrum) or 
with other geometric parameters for the loop. A thorough study of these coronal resonances 
for Alfven waves will be the subject of a future paper. 
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RTV scaling law 
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Fig. 1. — Temperature predicted by the RTV scaling law (see section I3T21) with respect to 
the obtained average temperature of the corona for the case of Alfven wave heating with 
monochromatic waves (panel a), different photospheric magnetic fields (panel b), different 
loop expansions (panel c), and different loop lengths (panel d). The values and corresponding 
symbols for each parameter are on the lower left side of each panel. 
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Nonlinear effects 

1.0L 




20 40 60 80 100 

Length along the loop [Mm] 

Fig. 2. — Nonlinear effects quantified through the ratio of the azimuthal velocity to the 
Alfven velocity, ({v^D / (v^) . Here we plot the average value in time of this quantity along 
the loop for different models. Panel a: base model, a loop of length 100 Mm, photospheric 
magnetic field of 2.3 kG and area expansion of 1000 between photosphere and corona. Panels 
b, c and d: Same as base model but with a photospheric magnetic field of lkG, with an area 
expansion of 500, with a length of 60 Mm, respectively. For each model three cases with 
different photospheric rms azimuthal velocity fields are plotted, the values of which are 
written in the upper right corner of each panel. The dotted lines correspond to a low value 
of (t> 2 ) 1//2 (leading to low coronal temperatures). The solid line corresponds to the value of 
(i;^) 1 / 2 for which the highest temperatures are found. The dashed line corresponds to a value 
for which the corona collapses during the simulation (see section [3.2. II) . 
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Monochromatic waves 
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Fig. 3. — Results of Alfven wave heating from monochromatic waves. The (longitudinal) 
waves have periods of 10, 25, 50, 100, 150, 200 and 300 s (corresponding symbols are written 
in the bottom left part of the figure). Average quantities are plotted with respect to the 
rms azimuthal velocity amplitude at the photosphere. Panels a, b d and e: temperature, 
density, plasma (3 and volumetric heating rate in the corona, respectively. The volumetric 
heating rate is calculated according to Eq. [121 Panel c: rms azimuthal magnetic energy in 
the photosphere. Panel g: maximum attained temperatures. Panels f and h: length and 
time span of the corona (with respect to total simulation time) , respectively. Blank symbols 
denote unstable coronae (see §3.2. lj) . 
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v s power map v. power map 
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Fig. 4. — Map of the power spectrum in logarithmic scale of the longitudinal velocity com- 
ponent v s (left) and azimuthal velocity component v<f> (right) for the first 20 Mm along a 
loop where 300 s monochromatic longitudinal waves are generated (from mode conversion of 
600 s Alfven waves). 
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Photospheric magnetic field 
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Fig. 5. — Results of Alfven wave heating with different photospheric magnetic fields. Fields 
of 2.3 kG (circles), 1.5 kG (triangles) and 1 kG (squares) are tested. See legend of Fig. [3]for 
explanation of the panels. 



-37- 



Loop expansion 
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Fig. 6. — Results of Alfven wave heating with different loop expansion factors A top /Ao, 
where A s is the cross-section area of the loop and s = 0, top correspond to the photosphere, 
apex of the loop (corona), respectively. Expansion of 1000 (circles), 600 (triangles) and 500 
(squares) are tested. See legend of Fig. [3] for explanation of the panels. 
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Fig. 7. — Results of Alfven wave heating with different loop lengths. Lengths of 100 Mm 
(circles), 80 Mm (triangles) and 60 Mm (squares) are tested. See legend of Fig. [3] for 
explanation of the panels. 
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Fig. 8. — Volumetric heating rate with respect to the considered parameters of the loop: 
period of the monochromatic longitudinal waves issuing from mode conversion of Alfven 
waves (top left panel), photospheric magnetic field (top right panel), expansion of the loop 
(bottom left panel), and loop length (bottom right panel). The volumetric heating rate for 
each case corresponds to its extremum values in the range between (^) 1/2 ~ 1.3 and 2.3 km 
s _1 found in the e panels of Figs. EJ El [6] and [71 



